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Приводятся результаты математического моделирования, методики и построения программ ЧПУ для обра- 
ботки деталей сложной геометрической формы в условиях ОАО «Роствертол». 
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Введение. Для авиационной промышленности характерно изготовление изделий, обладающих 
малой и изменяющейся жёсткостью по координатам перемещения инструмента относительно за- 
готовки, а также деталей сложной геометрии, при обработке которых меняется направление 
формообразующих движений. При обработке таких деталей приходится считаться с тем, что тра- 
ектории движения исполнительных элементов станка, программируемых в управляющей ЭВМ, 
отличаются от траекторий формообразующих движений инструмента относительно заготовки, по 
крайней мере, на величину упругих деформационных смещений вершины инструмента относи- 
тельно точки его закрепления в исполнительном элементе. Динамические структуры металлоре- 
жущих станков вместе с процессом обработки обладают следующими свойствами [1]: 

— Они являются диссипативными. Поэтому в пространстве состояния множеству И(Х) со- 


ответствует множество траекторий формообразующих движений, которое является притягиваю- 
щим, то есть соответствующие траектории формообразующих движений являются аттракторами. 

— Как правило, с учётом физической реализуемости, можно выбрать такой вектор управ- 
ления, формируемый в координатах исполнительных перемещений станка, которому соответству- 
ет асимптотически устойчивая траектория формообразующих движений, которая обеспечивает 
требуемые показатели геометрического качества обработки. 

Для согласования программы ЧПУ, траекторий исполнительных элементов и траекторий 
формообразующих движений инструмента относительно заготовки, в случае, если обеспечивают- 
ся требуемые показатели геометрического качества изделий, возможно два подхода. Первый, 
широко используемый в настоящее время, основан на принципах подчинения всех координат про- 
грамме ЧПУ. Второй основан на развиваемой в настоящее время синергетической теории управ- 
ления и на определении И(Х)) таким образом, что в управлении используются естественные ди- 


намические свойства системы, в том числе и связи, формируемой процессом резания. При этом 
траектории формообразующих движений, обеспечивающие необходимое геометрическое качест- 
во, формируются естественным образом в данной динамической структуре и являются аттракто- 
рами системы [1, 2]. Второй подход, основанный на использовании всех дополнительных связей, 
рассматриваемых в парадигме расширения пространства состояния, применяемой в синергетиче- 
ской концепции управления, позволяет фактически включить в управление все естественно суще- 
ствующие связи. 
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Пусть определено многообразие траекторий К, = 2, обеспечивающих требуемые показа- 
тели геометрического качества обработки. Ему соответствует многообразие траекторий исполни- 
тельных элементов №, е Х, и далее программ ЧПУ №, е И, то есть векторов управления И(Х). 


Если учитывать естественно существующие в динамической структуре связи и закономерности их 
изменения вдоль траектории движения исполнительных элементов, то именно на многообразии 
Х, должна строиться программа ЧПУ. Таким образом, для определения многообразия М, необ- 


ходимо решать обратные задачи динамики для уравнений, связывающих пространства Х, Хи Ц. 

Алгоритм определения траекторий формообразующих движений. В статье ограничимся 
случаем токарной обработки. Случаи фрезерования концевыми фрезами, рассверливания, обра- 
ботки отверстий многолезвийными инструментами подробно рассмотрены в монографиях [1—3]. 


Будем полагать заданными траектории исполнительных элементов, т. е. {Х, (Х,),аХ, (Х,)/ а}, 


/=1,2,3,4. Их можно представить в функции времени, т.е. {Х, (и (6}, где 


1 


И, (Е) =аХ, (Е) /аЕ . Каждой точке вектора Х соответствует пара векторов упругих деформаций ин- 


15 т 
струмента х = {х,,х,,х,} и заготовки у ={у,,У,,У.} . Координаты Х отсчитываются от центра 
вращения заготовки, находящегося на правом её конце (рис. 1). Упругие деформации 
ТР :Т. 
х= {ххх}, У={УУ,, Уз} отсчитываются по направлению, противоположному телу заго- 


товки (рис. 1). Таким образом, координаты формообразующих движений, определяющие геомет- 
рию детали, складываются из координат исполнительных элементов станка и упругих деформа- 


ций. Например, текущее значение радиуса г” в точке Х(”, определяющее значение диаметра 
Р®, равно 

г® (х’) = 2, (Х.)=Х, (ху) +Х, (ху) в (х’). (1) 
Тогда точка 2, (Х.) — координата формообразующих движений, а точка Х, (Х,} — точка коор- 


динат исполнительных элементов станка. Имеет также место упругое смещение точки контакта 


инструмента и заготовки по отношению к координате Х{” на величину х, (х’) +у, (х’). 


Выясним закономерности формирования упругих деформаций при заданных траекториях 


1 


{Х ,(Е),и (#)}. Упругие же деформации прежде всего зависят от сил. Силы резания при обработ- 


ке на металлорежущих станках формируются в результате пересечения тела режущего инстру- 
мента с заготовкой с учётом деформаций х, у (рис. 2). Уровень моделирования сил резания 


РЕ ЕЕ, в. в координатах состояния может быть различным. Для того чтобы выяснить ос- 


новные закономерности преобразования траекторий, будем анализировать медленные формооб- 
разующие движения, то есть полагать, что процесс резания по координатам упругих деформаций 
инструмента относительно заготовки является устойчивым и в пределах импульсной реакции под- 
системы инструмента и заготовки система является замороженной. Более точные обоснования 
иерархии систем дифференциальных уравнений динамики даны в монографии [1]. 
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Рис. 1. Система координат, определяющая состояние управляемой динамической системы точения без учёта электромеха- 
нического преобразователя 


Ось вращения шпинделя 





1,5 


Рис. 2. Схема формирования площади срезаемого слоя и координаты формообразующих движений по направлению Х, с 
учётом упругих деформаций 


Тогда при формировании сил резания А = а Е,, Е можно принять следующие гипоте- 


зы (рис. 2) [1]: 
— модуль сил пропорционален площади срезаемого слоя 5, т. е. 


Е =р(И, 5; (2) 
— справедлива гипотеза неизменной ориентации сил в пространстве, т. е. 
Е={Е, Е, Е} =Е {^^^}, (3) 


где Л = {^,, Л, к. — вектор угловых коэффициентов. Заметим, что коэффициент р, имеющий 


смысл давления на переднюю поверхность инструмента, также зависит от износа инструмента. 
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Таким образом, силы при неизменных условиях определяются величиной подачи на обо- 
рот $р (6) и глубиной резания & (#) (рис. 2) 


5(#)=$, (66 (6-60, (4) 

1 09099. : > ь 
где Е 55 (Е) = | [и (К)-и, (Е) -и,, (Е) — путь, пройденный вершиной инст- 

9ф< +9ф Е-т 
румента за время одного оборота Тс учётом скорости упругих деформаций в направлении 
Х} (ии, ). 

Заметим, что в (4) предполагается движение инструмента вдоль обрабатываемой детали. 
При его движении по направлению оси вращения заготовки необходимо учитывать срезание ма- 
териала на предыдущем обороте. Тогда на временном отрезке движений инструмента в сторону 
оси вращения детали величина глубины резания фактически определяется величиной подачи, то 


Е 
есть определяется по зависимости $, (Ё) = | [и (Ки, (Е -и, (Е) |. При движении инстру- 
ЕТ 


мента вдоль оси вращения заготовки (это наиболее типичный случай для токарной обработки) 
силы, формируемые процессом резания, определяются следующими интегральными преобразова- 
ниями: 


вр] | ео [5 ®-мож-ж 0-х) 


Е 2 (5) 
< (И -и,(0-и, а | 


т 


Причём знак скорости упругих деформаций у, И, зависит от знака ЧЕ, /4Ё . Заметим, 


Е 
что Г, (х,)- [м (К аЁ-х, (Е)- у, (Е)>0, в противном случае контакт инструмента с заготовкой 
0 


отсутствует и процесса резания нет. 
Введём в рассмотрение матрицы жёсткости подсистемы инструмента с“) Ч Г 


5,К=1,2,3 и заготовки с®) = [ку], 5,К =12,3. Во многих случаях бывает справедливым 


с® (Х, ) = сопз. Пренебрегая в (5) квадратичными формами относительно малых упругих де- 
формаций, а также учитывая, что © есть величина малая, получаем линеаризованное выражение 


для определения связи траекторий исполнительных элементов с траекториями упругих деформа- 
ций инструмента относительно заготовки: 
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Сл хХ1 + СолхХо + 


нь БаЕ-Ь | у, (6 - и 


Сл, хХ 
2х + С>2хХо + С, 2, хХ 
С: з,хХ1 + Созх Ха + С, 3, хХ 


= Л,р| &, Те Ба-Ь | у. (6% - | [% (си, даа 
- ЕТ С=Ё-Т &=0 
=Л,р ВХ 


(6) 
т Гб да 


ИЕ На, — Иа - | ПАСАТ а 


| 


Су зу У + Соз,уУ2 + Сзз,Уз = ль [и Иа - | Иа - | [и (о а 


- =Е-Т &=0 


а(у, +х,). а (у, +х:) 
где и, = Е И р 
В частности, если скорость И можно считать постоянной (усреднённой) в пределах оборо- 


та заготовки и величину припуска в пределах оборота также &, = соп$® ‚ то для вычисления упру- 





. В (6) полагается, что величина припуска &Ё, = соп5%. 


гих деформаций хи уна каждом обороте из (6) получаем 











сх" =Е, (7) 
где 
[сх +А,р5> ба» бах + М РЁ ^,р5> 0 Л,рЕ, 
РЭ Е: Зо: ОЕЕАОЬ Л›Р5> 0 Л.РрЕ» 
ря СКА рь С БЕ 2ОЕ Л-Р5> 0 Л-РЕ» 
Л,р$> 0 Л,рЕ, Су + Р5р С.4, Са, +МРЬЕ 
Л›р5> 0 Л.РрЕ» Силу +АР5р С.2у Сз2, +АРЁ 
|  ^Л-р5 0 Л.Р, Су + АЗР5р Созу Сзз + А3РЬ, | 
х@7) = Ук уз} } 
РА, [5 + ых, (Е-Т)+ у, (Е-Т) 
РА, [2>5> +Ёьх, (#-Т)+Ёьу.(Е-Т)| 
РА, [255> + ых, (#-Т)+ ву, (Е-Т) 
ыы РА, ть, $р +1х, (Е-Т)+ у, (#-Т) 
РА» [Е,5> +1 х, (Е -Т)+вьу, (Е-Т). 
) 











РА; | 6,5» +вх.(Е-Т)+Ьу, (#-Т)| 
Вектор КЕ кроме технологических режимов включает заданные функции х, [4 -Т) и 
у, [4 -Т) ‚ которые определяются на отрезке (0, #7) при определении деформаций в момент ЕЁ 


Если вычисление деформаций осуществляется в пространстве состояния по координате № от 
оборота заготовки к обороту, то х, (#-Т) и у, (#-Т) равны значениям упругих деформаций на 


предыдущем обороте. Таким образом, суммарная матрица жёсткости, имеющая размерность 
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бхб ‚, складывается из матриц жёсткости инструмента, заготовки и коэффициентов, зависящих от 
технологических параметров, а также параметров процесса резания (определяются обобщённым 
коэффициентом р). 

Кроме этого, будем считать направляющие станка идеальными, т. е. не вносящими иска- 
жений в траектории формообразующих движений. Принципиально кинематические возмущения 
со стороны направляющих можно учесть, если задана их статистика. Тогда для обеспечения, на- 
пример, неизменности диаметра необходимо выполнить условие 

Х, - (х, + у, ) = сое, (8) 
которое можно выполнить на основе решения обратной задачи динамики из (7) различным под- 
бором Х, (х,) и $, (х,) . Условие (8) фактически определяет один из инвариантов функциониро- 


вания станка в синергетической теории управления [4]. Кроме этого, ЕЁ, и $, в (7) есть парамет- 


ры управления. Подчеркнём, что траектории исполнительных элементов через изменения пара- 
метров управления смещают траекторию формообразующих движений. В общем случае пере- 
страивают динамическую систему резания. 

Связанная между собой совокупность Х, (Х.), 5, (Х,), изменяющаяся по координатам 
перемещения суппорта вдоль оси вращения заготовки, после элементарных преобразований даёт 
совокупность траекторий исполнительных элементов, обеспечивающих требуемое в данном слу- 
чае условие (8). Эти траектории и формируют многообразие №, е Х траекторий исполнительных 
элементов станка, обеспечивающее условие (8), т.е. неизменность диаметра — одна из целей 
процесса резания. 

Приведённые условия легко обобщаются на случай, когда деталь имеет сложную геомет- 
рическую форму. Здесь ограничимся случаем, когда геометрия детали допускает кусочно- 
линейную аппроксимацию (линейная интерполяция) и выделим на детали отдельный фрагмент 
(рис. 3). На этом фрагменте обработка ведётся с одновременным изменением скоростей продоль- 
ной и поперечной подачи (в данном случае обрабатывается элементарный конус). При анализе 


траекторий в пределах /го оборота частоту и, будем считать неизменной, т. е. 7® =1/щ). Если 
заданы скорости И/ (), И, (Е) и их вектор ориентирован по направлению образующей, то величи- 
на подачи на #м обороте будет равна 


Е (о 


ти 

Причём при линейной интерполяции для обеспечения движений по образующей М - № 

(рис. 3) должно выполняться требование И, (Е) (Е) = К“) = соп&Ё, если упругие деформации 

отсутствуют, то многообразие траекторий частот вращения двигателей исполнительных элемен- 

тов станка {и (Хы, (х,}}, обеспечивающих движение по желаемой образующей 
«М-М », определяется соотношениями: 

и, (Х,) = тг (Х4®); 


ый (10) 
м =" &,, 
где У, — скорость резания, &,, &, — частоты вращения двигателей поперечной и продольной 


подач. 
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Рис. 3. Схема формирования формообразующих движений при точении детали сложной геометрической формы 


Если И, постоянна и соответствует традиционно определяемой оптимальной скорости ре- 
зания, то частота вращения шпинделя есть функция радиуса, т. е. координаты Х, ; коэффициент 
км) определяется углом наклона образующей М -М (см. рис. 3). Упругие деформации инстру- 
мента в направлениях Х, и Х, искажают реальную образующую (рис. 4). Если упругие деформа- 


ции отсутствуют, то в процессе управления приводами за счёт асимметрии динамических харак- 
теристик приводов может нарушаться одно из условий (10), что не только вызывает отклонение 
реальной траектории относительно образующей М -— №, но и приводит к дополнительным возму- 
щениям. 

Для выяснения механизмов отклонения траекторий введём в рассмотрение вектор скоро- 


стей И“) = и (Жи (х,) ‚ отличающийся от И = М (Х,) ‚И, (х,}} матрицей перехода, 


1 3 


зависящей от угла ф“""), определяемого углом между осью вращения заготовки О,О' и обра- 

зующей ММ (рис. 4). Связь между скоростями определяется геометрическими соотношениями 
ИМ (Х)=ФИ(Х.), (11) 

И Фил а 

— 505 Фи, УПФил 


Если И, (Х, п фи, =И, (Х,) 60$ Фи, , то суппорт движется по прямой М-М№. Тогда фор- 


где о 


мирование сил резания соответствует ранее проанализированному случаю, но вместо скорости 
и, (Х.) необходимо рассматривать скорость ИГ" =и, (Х, ) япфии +И, (Х.)с05фи, и учитывать 
изменение угловых коэффициентов ориентации сил резания. Очевидно, что условию 
и, (Х, п фи =И (Х,)с0$ф„, соответствует выполнение второго соотношения в (10). Тогда для 


вычисления упругих деформаций для подсистемы инструмента (обрабатывается жёсткая деталь) 
можно воспользоваться соотношениями 
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сх Е, (12) 
С, 1 60$ Фил — С31 ЭПФии 65: С. ЯПФии + С; 1 60$ Фиу 
где с“) = с, со5фии — С: ЭП Фии С,, С, ЯПФии +6; > С05Фии |, причём коэффициенты с,,, 
С! з С05 Фин —С: 3 УПФии С.3 С; ЧПФии + С; ; 60$ Фи 
5,К =1,2,3 есть элементы матрицы жёсткости в системе координат х, т.е. при фи, =0 
с\) = с. Для определения вектора управления (программы ЧПУ) полученное многообразие 
траекторий №, необходимо перевести в многообразие №, . Эта операция осуществляется на 
основе прямой подстановки траекторий {Х (Е), аХ(Е )/ае} =Х, В уравнения управляемых дви- 


гателей исполнительных элементов. Тем самым учитывается не только особенность преобра- 
зования траекторий исполнительных элементов в траектории формообразующих движений, 
но и динамические особенности приводов, в частности, асимметрия их динамических харак- 
теристик. 


Заготовка 








Рис. 4. Связь площади срезаемого слоя с траекториями И, (Х,) , И, (Х,} и упругими деформациями х, (Х.), х, (Х,) 


Пример использования синергетической концепции. Рассмотрим модельный пример, соот- 
ветствующий точению вала диаметром а =20,0 мм и длиной [ = 430 мм из Ст45. На рис. 5 при- 
ведена расчётная схема, а на рис. 6 — диаграммы изменения элементов матрицы жёсткости под- 
системы заготовки, рассчитанные по методу конечных элементов. Матрицы жёсткости рассчитаны 
на основе программного комплекса ИСПА-386 — интегрированной системы прочностного расчёта 
методом конечных элементов. Для рассматриваемого модельного примера параметры матрицы 
жёсткости точек крепления заготовки и подвески режущего инструмента приведены в табл. 1. 
Для расчёта смещения точки равновесия х, + у, =и/ , усреднённой по периодам вращения заго- 
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товки, и определения траектории «медленного» изменения скорости подачи {Х, (т), аХ, И} ЕХ 


необходимо иметь значения параметра р, которые приведены в табл. 2 для случая обработки 
Ст45 трёхгранными неперетачиваемыми пластинами из твёрдого сплава Т15Кб при варьировании 


скорости резания в диапазоне И, = (60,0 - 180,0) м/мин для величины припуска ЕЁ, =2,0 мм. 





Рис. 5. Расчётная схема для определения матриц жёсткости заготовки по длине контакта инструмента и детали 






































Таблица 1 
Исходные матрицы жёсткости системы резания 
Наименование Са, С. >, С С Са, С. 2, 
кг/мм кг/мм кг/мм кг/мм кг/мм кг/мм 

Матрица жёсткости крепления | 8,0.10“ 6,8.10* 5,8.10* 4,2.103 2,8.103 2,4.103 
заготовки в патроне [6] 
Матрица жёсткости крепления | 13.10“ 0,8.10* 7,2.10* 0,5.103 2,8.103 2,3.103 
заготовки в задней бабке 
[< 
Матрица жёсткости подсисте- | 0,5.10“ 0,3.10* 0,3.10* 0,08.10* 0,05.10* 0,04.10* 
мы инструмента [6% | 

Таблица 2 


Изменение коэффициента р при варьировании скорости резания 























У», м/мин | р, кГ/мм? 
60,0 290,0 
100,0 250,0 
140,0 220,0 
180,0 210,0 








На рис. 7 приведены рассчитанные траектории изменения скорости подачи при различных 
значениях скорости резания У,, при которых выполняется условие и/ = сопзЁ. Во всех случаях 


величина припуска считается постоянной равной #, = 2,0 мм. На рисунке вместо скорости приве- 


дена частота вращения шпинделя, которая, что очевидно, при неизменном диаметре детали от- 
личается от скорости на коэффициент. 
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6) 
Рис. 6. Изменение матриц жёсткости заготовки вдоль координат обработки: 
а — диагональные элементы; 6 — недиагональные элементы 


Все траектории, принадлежащие поверхности, приведённой на рис. 7, и удовлетворяющие 
условию @Х, /4ё >0 , обеспечивают обработку с постоянными значениями упругих деформаций 


инструмента относительно заготовки в направлении, нормальном к оси вращения шпинделя. На 
приведённой иллюстрации траектория А соответствует обработке при неизменной частоте вра- 
щения шпинделя, а траектория В — при монотонно возрастающей. Направление движения суп- 
порта показано стрелками. На участке 1—2 скорость остаётся неизменной, она соответствует ско- 
рости врезания инструмента в заготовку. Кроме этого, по расчётам на участке врезания наблюда- 
ется резкое возрастание скорости подачи, которое на иллюстрации не показано. Оно длится в 
течение 2—3 оборотов и связано с переходными процессами в зоне резания. Сразу обратим вни- 
мание на многовариантность таких траекторий. 

Естественно поставить вопрос о выборе на многообразии некоторой единственной траек- 
тории, являющейся оптимальной в определенном смысле. Методика выбора траектории опти- 
мальной скорости резания вдоль траектории движения инструмента относительно заготовки по 
критерию минимизации приведённых затрат на изготовление партии изделий изложена в работах 
[5, 6]. Здесь же подчеркнём, что речь идёт не об оптимальной, но неизменной скорости резания, 
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традиционно рассматриваемой в технологических исследованиях, а 0б оптимальной траектории. 
Такая траектория вычисляется на основе решения неклассических вариационных задач, в данном 
случае — на основе принципа динамического программирования. 
6, , об/мин 
3000 


2000 


1000 


2,5 





И, „мм/с 


Рис. 7. Многообразие траекторий исполнительных элементов станка, обеспечивающее постоянство упругой деформации 
инструмента относительно заготовки при изменениях скорости подачи и частоты вращения шпинделя 


Заключение. Опыт использования синергетической концепции при обработке изделий в услови- 
ях ОАО «Роствертол» показал возможность снижения трудоёмкости изготовления изделий ответ- 
ственного назначения на 30—40 %. Однако внедрение этого метода в полной мере требует мо- 
дернизации существующих на предприятии систем ЧПУ. Сущность модернизации заключается в 
установке сопроцессоров непрерывного управления траекториями исполнительных элементов, 
позволяющих обеспечивать изменение в каждом кадре управления скоростью по перемещению. 
При построении программы ЧПУ используется линейная интерполяция, но она осуществляется не 
по геометрическому образу детали, а по многообразию №, . При этом используются связи, естест- 


венно существующие в динамической системе станка и формируемые в процессе резания, кото- 
рые автоматически обеспечивают принцип когерентности при управлении. Изложенный подход 
можно распространить на обработку концевыми фрезами, а также на другие процессы обработки, 
в том числе многолезвийными инструментами [1]. Этот подход, дополненный сопроцессорами ди- 
намического мониторинга обработки [7—9], характеризует новый этап в создании систем управ- 
ления станками в интегрированном производстве. 
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